perm filename TABLE.SAI[1,VDS] blob
sn#498112 filedate 1980-02-03 generic text, type C, neo UTF8
COMMENT ⊗ VALID 00002 PAGES
C REC PAGE DESCRIPTION
C00001 00001
C00002 00002 BEGIN "TABLE"
C00006 ENDMK
C⊗;
BEGIN "TABLE"
REAL ARRAY A[1:3,1:6],B[1:3,1:6],R[1:3],T[1:3,1:3],L[1:3,1:6],S[1:6];
REAL LI,A1,B1,G1,X1,Y1,Z1,DX,DY,DZ,ALTHA,BETA,GAMMA,PI;
DEFINE CRLF="'12&'15";
PI←3.14159265
OUTSTR ("TYPE IN A1,B1,G1 ");
INCHWL A1,B1,G1;
READ (5,205)X1,Y1,Z1
WRITE (6,220)
WRITE (6,230) A1,B1,G1
WRITE (6,235)
WRITE (6,240) X1,Y1,Z1
DO 5 J=1,6
READ (5,205) (A(I,J),I=1,3)
5 CONTINUE
DO 10 J=1,6
READ (5,205) (B(I,J),I=1,3)
10 CONTINUE
ALPHA=0
BETA=0
GAMMA=0
READ (5,205) DX,DY,DZ
READ (5,200) ALPHA,BETA,GAMMA
WRITE (6,245)
WRITE (6,250) DX,DY,DZ
WRITE (6,255)
WRITE (6,260) ALPHA, BETA,GAMMA
R(1)=DX+X1
R(2)=DY+Y1
R(3)=DZ+Z1
ALPHA=(ALPHA+A1)*(2.*PI)/360.
BETA=(BETA+B1)*(2.*PI)/360.
GAMMA=(GAMMA+G1)*(2.*PI)/360.
C1=COS(ALPHA)
C2=COS(BETA)
C3=COS(GAMMA)
S1=SIN(ALPHA)
S2=SIN(BETA)
S3=SIN(GAMMA)
T(1,1)=C3*C2
T(1,2)=(C3*S2*S1)-(S3*C1)
T(1,3)=(C3*S2*C1)+(S3*S1)
T(2,1)=S3*C2
T(2,2)=(S3*S2*S1)+(C3*C1)
T(2,3)=(S3*S2*C1)-(C3*S1)
T(3,1)=-S2
T(3,2)=C2*S1
T(3,3)=C2*C1
DO 60 J=1,6
DO 50 I=1,3
L(I,J)=(T(I,1)*A(1,J)+T(I,2)*A(2,J)+T(I,3)*A(3,J))+R(I)-B(I,J)
50 CONTINUE
60 CONTINUE
DO 70 J=1,6
S(J)=SQRT(L(1,J)*L(1,J)+L(2,J)*L(2,J)+L(3,J)*L(3,J))
70 CONTINUE
WRITE (6,265)
WRITE (6,270)
DO 80 J=1,6
WRITE (6,275) J,S(J)
80 CONTINUE
200 FORMAT (3F6.3)
205 FORMAT (3F5.3)
220 FORMAT (' ','THE COMPONENTS OF INITIAL ANG. POSITION IN DEG. ARE')
230 FORMAT (' ','ALPHA=',1F7.3,2X'BETA=',1F7.3,2X,'GAMMA=',1F7.3)
235 FORMAT ('0','COMPS OF THE INITIAL POS. OF MOV FRAME ORIGIN ARE')
240 FORMAT (' ','X=',1F6.3,2X,'Y=',1F6.3,2X,'Z=',1F6.3)
245 FORMAT ('0','FOR DISPLACEMENTS FROM THE INITIAL POS. OF;')
250 FORMAT (' ','X=',1F6.3,2X,'Y=',1F6.3,2X,'Z=',1F6.3)
255 FORMAT ('0','AND COMPONENT VECTORS OF;')
260 FORMAT (' ','ALPHA=',1F7.3,2X,'BETA=',1F7.3,2X,'GAMMA=',1F7.3)
265 FORMAT ('0','THE LENGHTS OF THE STRAIN ELEMENTS WILL BE;')
270 FORMAT ('0',T30,'ELEMENT NUMBER',10X,'LENGTH,IN.')
275 FORMAT (' ',T38,I1,18X,1F6.3)
STOP
END